l/f a noise from correlations between avalanches in self-organized criticality 
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We show that large, slowly driven systems can evolve to a self-organized critical state where long 
range temporal correlations between bursts or avalanches produce low frequency 1/ f a noise. The 
avalanches can occur instantaneously in the external time scale of the slow drive, and their event 
statistics are described by power law distributions. A specific example of this behavior is provided 
by numerical simulations of a deterministic "sandpile" model. 
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The ubiquity of 1/ f a noise in nature is one of the 
oldest problems in contemporary physics still lacking a 
generally accepted explanation, despite much effort. The 
phenomenon is characterized by a decay with a non- 
trivial awl found in the power spectrum S(f) of a given 
time signal at low frequencies /. It has been observed in 
a huge number of diverse systems, many of which are far 
from equilibrium and show an avalanche- like dynamics. 
Examples include earthquakes ffijM, combustion fronts 
[@], chemical reactions Q, flux motion in superconduc- 
tors H , and Barkhausen noise || , to name only a few. 

In the search for a general dynamical mechanism for 
1// Q noise in far from equilibrium systems, Bak, Tang 
and Wiesenfeld (BTW) proposed the concept of self- 
organized criticality (SOC) [0j8|. This refers to the ten- 
dency of spatially extended, slowly driven systems to 
organize into a state with fractal spatial and temporal 
properties that is also characterized by self-similar dis- 
tributions of event (avalanche) sizes. Although the orig- 
inal "sandpile" model, which was introduced to exem- 
plify the SOC concept, exhibits a scale-free distribution 
of avalanches, its noise spectrum is of the form 1/f 2 [[|. 
Variants of the original sandpile model do indeed ex- 
hibit nontrivial l/f a noise, but in these cases jl(J the 
avalanche event distributions are not critical, implying, 
possibly, mutual exclusivity between a SOC mechanism 
and mechanisms giving long range temporal correlations 
such as 1/ f a noise. 

In fact, it has recently been argued in the context of 
solar flares, transport dynamics in magnetic confinement 
devices, and other areas, that the presence of tempo- 
ral correlations between events excludes SOC as an un- 
derlying mechanism. For instance, in Ref. pd| ] it was 
stated that for SOC "one expects no correlation between 
... bursts". See Ref. [Oj for a discussion and references. 
The argument entails a narrowing of SOC to the phe- 
nomenology of certain, specific "sandpile" models, which 
is in sharp contrast to the original idea (?]||] , and is erro- 
neous as we show here. 

In order to demonstrate explicitly that SOC can pro- 
vide a dynamical mechanism which gives correlations be- 
tween bursts leading to 1/ f a noise, we study a slowly 



driven, deterministic "sandpile" model. It exhibits a 
power-law distribution of avalanches, as well as l/f a 
fluctuations in the pattern of dissipation over the slow 
temporal domain of the external drive. Despite the fact 
that the time scales of the driving and of the avalanche 
events are completely separated, the critical behavior of 
the power spectrum is solely determined by the critical 
properties of the avalanche size distribution. In addition, 
the results we find are robust with respect to changes in 
the definition of the time scale associated with the driv- 
ing. These observations constitute a proof that a SOC 
mechanism can give low frequency 1//" noise due to 
correlations between power-law distributed avalanches, 
without imposing temporal correlations in the external 
drive. 

Actually several different models of SOC, describing 
e.g. traffic [|l3| and evolution J|,[l4| do exhibit nontriv- 
ial l/f a noise over the temporal domain of individual 
avalanches. In this context, it is important to distinguish 
between temporal correlations within the time span of 
individual avalanche events, and low frequency noise ob- 
served over a much longer temporal regime of an arbitrar- 
ily slow external drive. Measured in this external time 
scale, the individual avalanches can occur almost instan- 
taneously. This is a situation often encountered in Na- 
ture. In this case 1/ f a noise must arise from correlations 
between avalanches. These correlations can either be in- 
duced by long range temporal correlations in the external 
drive, or be an intrinsic part of the self-organization pro- 
cess itself. Sanchez et al have found long range temporal 
correlations in SOC when the "sandpile" is driven by an 
external source of low frequency noise |T^ ] . 

However, it is important to clarify if SOC it- 
self can spontaneously generate both critical avalanche 
statistics and long-range temporal correlations between 
avalanches, in the presence of a temporally uniform, slow 
external drive. Surprisingly, this question has not yet 
been decided, despite the fact that, as in the case of 
earthquakes, a scale-free avalanche event dynamics (e.g. 
the Gutenberg-Richtcr law) can be observed in many far 
from equilibrium phenomena which exhibit l/f a noise 
(and non-Poissonian inter-event statistics |l|,^5|) at vastly 
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longer time scales (e.g. days to years) compared to the 
individual events (e.g. seconds). Perhaps the lack of a 
sufficiently convincing "sandpile" model has contributed 
to the ongoing debate about the relevance of SOC to so- 
lar flares, earthquakes, transport in tokamaks, and many 
other burst-like phenomena. 

The "sandpile" model we discuss was introduced by 
de Sousa Vieira to describe avalanches in stick-slip phe- 
nomena . It is close to the original array of connected 
pendula first discussed by Bak, Tang, and Wiesenfeld. In 
detail, the model is defined as follows: Consider a one- 
dimensional system of size L where a continuous vari- 
able f n > is associated with each site n (representing 
the force on that site). Initially, all /„ have the same 
value which is below a threshold fth- The basic time 
step of the driving force consists in changing the value of 
the force on the first site according to f\ ~ f t h + Sf 
with a fixed Sf. This can be considered as a slow 
external driving and leads to a fast relaxation process 
(avalanche) within the system. This relaxation consists 
of a conservative redistribution of the force at sites with 
fn > fth (toppling sites) according to f n = $(/„ - f th ) 
and f n±1 = f n±1 + A/ n /2. Here, A/„ is the change of 
force at the overcritical site and <!> a periodic, nonlinear 
function. The relaxation continues until all sites are sta- 
ble again, i.e., f n < fth for all n. Then, the driving at the 
first site sets in again. This definition is complemented 
by open boundary conditions, i.e., force is lost at both 
boundaries. Without loss of generality, we use a sequen- 
tial update and set f t h = 1, Sf = 0.1, = 1 — a[x] 
where [x] denotes x modulo 1/a, i.e., a sawtooth function, 
and a = 4. It was shown that the model evolves into a 
state of SOC where the avalanche distributions are scale 
free, limited only by the overall system size [Q. Note 
that like other sandpile models, the toppling rules are 
conservative, and the total amount of force in the system 
can only be changed at the boundaries. 
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FIG. 1. Fluctuations of the total force in the self-organized 
critical state for a system of size L = 1024. 



An appropriate choice of a time signal to detect 1/ f a 
noise, which takes the time scale separation between 
the slow external driving and the individual avalanches 
into account, is the total force in the system after each 
avalanche: 

L 

X(i) = J2fnd), (1) 

71=1 

where i is the avalanche number. This signal is shown in 
Fig. [I] and directly reveals the stick-slip character of the 
dynamics. 

Analyzing the power spectrum of X(i), we find a clear 
1/ f a decay with a cut off at low frequencies that shifts 
to even lower frequencies with increasing system size (see 
Fig. |^). In particular, a data collapse, shown in the inset, 
reveals the following scaling behavior: 

s(f) - pK^), (2) 

with a = 1.38 and (3 — 1.2. Here, g is a scaling func- 
tion that is constant for small arguments and decays as 
x~ a for large arguments. These results are somewhat 
dependent on the details of the definition of the model 
]l7t , a fact which does not change our main conclusion 
regarding SOC as a mechanism leading to 1/ f a noise. 
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FIG. 2. Power spectrum of the fluctuations in the total 
force for different system sizes. Note that the cut off moves to 
lower frequencies for larger system sizes. Inset: Data collapse 
with a = 1.38 and (3 = 1.2. 

To understand how the long range temporal correla- 
tions arise in this model, we have analyzed the fluc- 
tuations in the force at a single site, i.e. X n (i) — 
f n (i). It turns out that the local power spectra are sim- 
ple Lorentzians with a characteristic frequency that de- 
creases with increasing distance from the driving site. 
Thus, different time signals of the same system, such as 
the local and the total force signal, behave in a totally dif- 
ferent manner, confirming earlier results J18|. Moreover, 
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the observation of Lorentzians suggests that the l/f a 
noise could occur due to a superposition of independent 
and purely random processes each described by a char- 
acteristic time scale t c (see, e.g., Ref. However, the 
fluctuations of the force at different sites are correlated. 
This can already be deduced from the fact that the force 
at a given site can only change if all sites closer to the 
first site have discharged as well. A more careful analysis 
shows that the assumption of uncorrelated (local) signals 
would give a value of a far away from the observed value 
p9[ . Consequently, the long range temporal correlations, 
implied by the occurrence of 1/ f a noise, are stored in the 
spatial correlations embedded in the whole system. 

This is further confirmed by the fact that scaling re- 
lations connect a and (3 to the critical exponent charac- 
terizing the avalanche distributions. The quantity L 13 in 
Eq. (2) describes the scaling of the temporal cutoff and, 
hence, the scaling of the number of avalanches before 
events become uncorrelated. Thus, it can be related to 
the avalanche distribution. The distribution of avalanche 
sizes (the number of toppling events in the avalanche) is 
distributed as P(s) ~ s~ T G(s/ L D ), where r ~ 1.54 is the 
so-called histogram exponent, and the avalanche dimen- 
sion D ~ 2.20 gives the cutoff for the largest avalanche 
in a system of size L Jl6| ]. Avalanches in the power-law 
regime do not extend through the entire system; therefore 
some sites do not topple and the system retains memory 
of its previous force. However, avalanches larger than the 
cutoff, which are the ones that entirely span the chain, 
decorrelate the system because all sites topple and get 
a new force. The frequency of the large system wide 
events that decorrelate the force in the system scales as 
L -D(r-i) t]iuSj p _ £>( r _ i). Also, a scaling relation 

(s) ~ L was found numerically in |l6| ], and is a general 
result for many boundary driven SOC systems J20|, giv- 
ing the result that D(2 — t) = 1. Combining the two 
equations, (3 = D — 1 = (t — l)/(2 — r), which agrees 
very well with numerical results. 

Using scaling relations together with some previously 
obtained results concerning universality classes, we can 
relate the exponent a to the exponent r or D via the 
variance in the total force in the system: 



J S(f)df = v 2 =/(j2 M) - (E f»® 




(3) 



where the integral is over all frequencies. Since the expo- 
nent a > 1, the power in the signal diverges as L —> oo 
as L^( Q_1 ). As argued and supported numerically in Ref. 
JT(| , the "sandpile" model is in the same universality class 
as the original Burridgc-Knopoff train model. This lat- 
ter model has been conjectured to be in the universality 
class of interface depinning with the interface pulled at 
one end pp| - a result which is also supported by numer- 
ical simulations. Thus, the fluctuations in the force in 
the system have the dimension of fluctuations of force in 



the interface depinning problem, e.g. (J Q dxV 2 H(x)) 2 
(see Ref. |2C|] for more explanation). Since in that case, 
the height of the interface H has the dimension of L , 
then cr 2 ~ L 2 ( D ~ 2 \ Using the previous relation for /3, we 
get a = (3D-5)/(D-l) = (5t-7)/(t- 1). This again 
agrees very well with our numerical simulation results 
and is consistent with the results in ]2~lj] . The existence 
of these scaling relations makes it clear without ambigu- 
ity that, in this model, the critical avalanche dynamics 
and the long range temporal correlations belong insepa- 
rably together. 

Instead of using the avalanche number, we can choose 
different definitions for the slow time scale of the model. 
Since the force on the first site is set to fth + 5f to induce 
a new avalanche, one can choose to identify the amount 
of force added to the first site with the temporal inter- 
val between avalanches. This corresponds to a uniform 
driving. Such a choice leads to non-uniform time inter- 
vals between successive relaxation events which could in 
principle alter the behavior. The time signal defined in 
Eq. (|l|) is modified and becomes a sawtooth function 



x(t) = ]T /„(*) - 1 - u + £ uu), 



(4) 



where U < t denotes the time of occurrence of the last 
avalanche. The power spectrum of this time signal shows 
exactly the same l/f a decay as the signal defined in Eq. 
(jl]) as can be deduced from Fig. 0. There, we show 
the power spectrum of the differentiated signal X'(t) - 
called a dissipation signal for obvious reasons — which 
follows a power law with S(f) oc f 2 ~ a . In particular, 
the cut off to white noise behavior (a = 0) for very low 
frequencies occurs almost at the same frequency as before 
and the scaling with L is also unchanged. Note that 
white noise, (S(f) ~ /°) in X(t) in Fig. 2, corresponds 
to S(f) cx f 2 in X'{t) in Fig. 3. 

The dissipation signal resembles the form of a pulse 
train, i.e., 



X'{t) = J2hiS{t-U), 



(5) 



where hi is the dissipation in the total force due to the 
ith avalanche. This form of the signal is especially well- 
suited to clarify the source of the long time correlations. 
For instance, substituting hi by a constant should re- 
veal the correlations induced solely by the fluctuations in 
the time intervals between dissipation events. As shown 
in |l 22 1, correlations between subsequent time intervals 
can, indeed, lead to 1//" noise. However, here we find 
that this "return signal" is uncorrelated (see Fig. ||) ■ The 
absence of correlations in the time intervals together with 
the existence of a finite second moment of their distribu- 
tion due to the boundedness of the intervals by f t h + Sf 
implies that the power spectrum is the same as before 



3 



for frequencies below a certain, fixed fh ■ Hence, the long 
time scales are not affected by the fluctuations in the time 
intervals. In particular, this is true for any distribution 
of time intervals with finite second moment if the time 
intervals are independent of one another (central limit 
theorem). Thus, the l/f Q noise is due solely to the fluc- 
tuations in hi. We conclude that the long range temporal 
correlations in the model arc exclusively encoded in the 
state of the whole system and are not destroyed by the 
fluctuations due to changes in the definition of the time 
scale of the external driving. 



demonstrate that a SOC mechanism, with a power-law 
distribution of avalanches, can lead to long time corre- 
lations between avalanches and thus low frequency 1/ f a 
noise in slowly driven systems. 

We thank P. Bak for useful discussions. J. D. would 
like to thank the DAAD for financial support and the 
Imperial College for support and hospitality. 
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FIG. 3. Power spectrum of the different, differentiated, 
time signals for L = 256 (see text) . Note that in the avalanche 
and dissipation signal, an increase with f 2 as observed for the 
lowest frequencies corresponds to white noise in the integrated 
signal. The avalanche signal is shifted down by 8 decades. 

The fact that the occurrence of l/f a noise in the 
model cannot be attributed to a low-dimensional dynam- 
ics shows that the SOC mechanism is totally different 
and can be well-distinguished from other deterministic 
mechanisms which attempt to explain flicker noise as a 
chaotic phenomenon (see, e.g., Refs. psfl). Unlike those 
mechanisms, the dynamics discussed here cannot be re- 
duced to a renewal process with a power-law distribution 
of waiting times or step sizes. 

Finally, the temporal variations in the sizes of the 
avalanches show only trivial long time correlations. Fig. 
|3| shows the power spectrum of the avalanche signal, de- 
fined in (||) where hi is now the number of topplings in the 
ith avalanche. Clearly, a flat spectrum can be observed 
with a change to f 2 behavior at low frequencies due to the 
finite system size. In terms of the integrated signal X(t), 
this corresponds to a trivial 1/f 2 decay crossing over to 
white noise (/°) behavior at low frequencies. Hence, the 
underlying process can be considered as a random walk 
confined by the system size. 

To summarize, the interplay between building up of 
force from the external driving when the avalanches are 
small, balanced by dissipation of force at the bound- 
aries when the avalanches are larger is responsible for 
the flicker noise in this SOC model. Our results explicitly 
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